First-principles studies of water adsorption on graphene: The role of the substrate 
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We investigate the electronic properties of graphene upon water adsorption and study the influ- 
ence of the Si02 substrate in this context using density functional calculations. Perfect suspended 
graphene is rather insensitive to H2O adsorbates, as doping requires highly oriented H2O clusters. 
For graphene on a defective Si02 substrate, we find a strongly different behavior: H2O adsorbates 
can shift the substrate's impurity bands and change their hybridization with the graphene bands. 
In this way, H2O can lead to doping of graphene for much lower adsorbate concentrations than for 
free hanged graphene. The effect depends strongly on the microscopic substrate properties. 



Graphene, i.e. a monolayer of graphite, is the first 
truly two dimensional (2D) material [r, ^ and a promis- 
ing candidate for silicon replacement in semiconductor 
industry |3( or gas sensing applications [J, Since 
graphene's discovery the water adsorbates have been dis- 
cussed as impurities leading to doping [1,4], while chang- 
ing the electron mobility in graphene only surprisingly 
little. Up to now, the microscopic mechanism of this dop- 
ing without significant changes in electron mobility has 
remained unclear. Density functional theory calculations 
on single water molecule adsorbates on perfect free stand- 
ing graphene were in line with previous studies on 
carbon nanotubes (CNT) 0] and found H2O physisorp- 
tion but no H2O induced impurity states close to the 
Fermi level. Therefore, the doping effects found experi- 
mentally y are very likely due to more complicated 
mechanisms than interactions of graphene with single wa- 
ter molecules. The experiments dealing with the effect of 
water on graphene were carried out using graphene on top 
of Si02 substrates. In addition, for finite concentrations 
of H2O on graphene H2O clusters might form. Despite 
its importance for understanding the doping experiments 
and for possible graphene based applications, a mi- 
croscopic theory on the role of the substrate and H2O- 
clusters in the H2 0-graphene-interplay has been lacking, 
so far. 

In this letter we study the substrate and cluster for- 
mation effects in the water-graphene-system by means 
of density functional theory (DFT). We show that both, 
highly oriented water clusters as well as water adsorbates 
in combination with a defective Si02 substrate can lead 
to doping of graphene. To this end, we consider model 
systems (see Fig. [1]) for water and ice in different con- 
centrations of free standing graphene as well as for water 
interacting with defective Si02 substrates. By analyzing 
the involved dipole moments and comparison to electro- 
static force microscopy [s], we show that the Si02 sub- 
strate is crucial for obtaining doping by H2O adsorbates 
on graphene. 

In our DFT calculations, the electronic and structural 
properties of the graphene-substrate-adsorbate systems 



are obtained using generalized gradient approximation 
(GGA)[i,[l^ to the exchange correlation potential. For 
solving the resulting Kohn-Sham equation we used the 
Vienna Ab Initio Simulation Package fVASP) [TH with 
the projector augmented wave (PAW) [ill, [13] basis sets 
implemented, there. The graphene-substrate-adsorbate 
systems are modeled using supercells containing up to 83 
substrate atoms (Si, O and H), 12 adsorbate atoms and 
32 C atoms. 

Firstly, water adsorption on free standing graphene 
with different water concentrations is considered. To 
model a single H2O molecule on graphene, 3x3 graphene 
supercells were used. Full relaxation of H2O with oxy- 
gen or hydrogen nearest to the graphene yielded ad- 
sorbed configurations with binding energies of 40meV 
and 36meV and molecule sheet distances of 3.50 A and 
3.25 A, respectively. These values are in the same order 
as those obtained by Leenaerts et. al. indicating ph- 
ysisorption of single water molecules on graphene. Ana- 
lyzing the density of states/band structures for these two 
adsorption geometries we find qualitative agreement with 
Ref. [6] : None of these configurations exhibits energy lev- 
els due to the adsorbate near the Dirac point, as shown 
in Fig. [2] a) for the configuration with oxygen closest to 
graphene. The HOMO of H2O is more than 2.4 eV below 
the Fermi energy and its LUMO more than 3 eV above it. 
The absence of any additional impurity level close to the 
Dirac point shows that single water molecules on perfect 
free standing graphene sheets do not cause any doping. 

The supercell applied here, corresponds to an adsor- 
bate concentration of n = 2nm~^, which is well inside 
the range of concentrations (1-10 nm~^) found experi- 
mentally in Ref. The lateral dimension a = 4.5 A of 
the hexagonal ice Ih (OOOl)-surface unit cell corresponds 
to a concentration of 5.7 H20nm~^ per layer. Thus, in- 
creasing the H2O concentration significantly above the 
n = 2 nm~^ from above leads to water clusters or ice like 
structures, rather than isolated molecules. 

To gain insight into doping of graphene by water clus- 
ters and ice overlayers we studied fully relaxed bi-layers 
and four layers of ice Ih adsorbed on graphene. These 
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Figure 1: (Color online) Model systems for water interacting 
with graphene. a-c: Free standing graphene with single water 
adsorbates (a), a bilayer (b) and a tetralayer (c) of ice Ih. 
Carbon atoms yellow, oxygen red and hydrogen small blue 
balls, d-i: graphene on top of Si02 with every second (d-f) 
or eighth (g-i) surface Si atom forming a Q§ defect. Water 
adsorbates are considered on top of graphene (e) , and between 
graphene and the substrate (f), (h), and (i). Fully coordinated 
Si atoms are depicted as big blue balls, Si atoms at Q3 defects 
in green. 

overlayer structures have been proposed as basis of ice 
growth on various hexagonal metal surfaces [lH and 
can be modeled as {VSx y/S) R30° overlayer on the simple 
graphene unit cell. The lattice mismatch in this configu- 
ration is 0.23 A — on the same order as found for water 
overlayers on Ni(lll) [14] — and therefore a reasonable 
starting point for studying ice on graphene. 

The supercell bandstructures (Fig. [2]b) and c)) show 
that the electric field by proton-ordered ice on top of 
graphene changes the energy of graphene's nearly free 
electron bands. In contrast to pristine graphene, where 
these bands start at 3eV above the Dirac point or in 
the case of single H2O adsorbates on graphene (Fig. [2] 
a)), their bottom is at 0.6eV above and O.leV below the 
Dirac point for a bi- and a tetralayer of ice Ih on top of 
graphene. This shift is due to electrostatic fields induced 
by the II2O dipole moments and results in hole doping 
for the tetralayer of ice on graphene. 

The water adlayers cause a change in contact potential 
S(j)^ which can be estimated using the H2O dipole mo- 



ment p = 6.2 X 10~^^ Cm and the relaxed structures to 
be 5<p = 1.4V and 5.4V for a water bi- and tetralayer, re- 
spectively. While only the latter structure causes doping, 
the corresponding change in contact potential exceeds the 
experimental value, Pexp = 1.3 V, from Ref. Q by more 
than a factor of 4. However, water strongly diluted in 
N2 has been found to cause hole doping in graphene on 
Si02[4]. Given these two experiments, doping due to 
multiple fully oriented ice overlayers as in Fig. [2] is likely 
not the most important interaction mechanism for water 
and graphene. 

We now turn to studying the effect of the Si02 sub- 
strate in the water-graphene interplay. The experiments 
with graphene on top of Si02 used substrates, which were 
created by plasma oxidation of Si. The Si02 surface 

created in this way is amorphous and its electronic, struc- 
tural and chemical properties are challenging to model 
from first principles. To obtain, qualitative insight to 
the most important physical mechanisms it is, however, 
a reasonable starting point to consider crystalline Si02 
in the /3-cristobalite form as a substrate [16]. 

The (111) surface of this modification can be used to 
create the most likely defects on Si02 amorphous sur- 
faces: These are so called Qs and Q\ defects [l^ hav- 
ing one under coordinated silicon and oxygen atom, re- 
spectively. Furthermore, the unit cell of this surface is 
nearly commensurate with the_graphene lattice: The lat- 
tice constant asi02 = 7.13 A [18] results in a surface unit 
cell asi02/V^ = 5.04 A, which is 4% larger than twice 
the lattice constant 2ao = 4.93 A of graphene. As the 
Si04 tetrahedra in Si02 are known to adjust to exter- 
nal strain easily, we model graphene on Si02 by 2 x 2 
or 4 X 4 graphene supercells with lateral dimension 2ao 
or 4ao. As substrate we put 4% laterally strained and 
hydrogen passivated /3-cristobalite with 6 Si atoms per 
surface unit cell in vertical direction. We then created 
the defects by removing H passivation atoms, added the 
H2O adsorbates and relaxed until all forces were less than 
0.08eV-A-i. In this way, we consider passivated and de- 
fective Si02 surfaces — the latter containing either un- 
dercoordinated silicon or oxygen atoms — as substrate 
for graphene. The effect of water exposure is simulated 
by putting water molecules on top of graphene as well as 
between graphene and the substrate. 

Graphene on top of fully passivated Si02 means two 
inert systems in contact with each other. Consequentely, 
there are no bands in addition to graphene's Dirac bands 
at the Fermi level and no doping occurs. (The band 
structure is not shown here, for brevity.) This changes 
strongly for graphene on defective Si02. As a model sys- 
tem, we study defects in /3-cristobalite (lll)-surfaces. 
Depending on the supercell size, 2x2 and 4 x 4, in these 
periodic structures every second and eighth surface Si 
atom is under coordinated, respectively, and forms a 
defect. (See Fig. [1] (d-i).) 

These defects lead to additional states in the vicinity 
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Figure 2: (Color online) Bandstructures of supercells with fully relaxed single molecules (a), a bilayer (b) and a tetralayer (c) 
of water on graphene, corresponding to Fig. [T]a-c), respectively, are shown. The graphene tt bands are marked in green, the 
nearly free electron bands in blue. Due to the H2O dipole moments, graphene's nearly free electron band is shifted with respect 
to its TT bands. 



(±leV) around the Fermi level. (See Fig. [3] a) and d).) 
The avoided crossing in Fig. [3] a) indicates significant 
hybridization of the defect and the graphene bands and 
demonstrate the impurity state's importance for conduc- 
tion electron scattering in graphene. In the pure Si02 
graphene model systems (Fig. [3] a) and d)), the impurity 
bands do not cross the Fermi level. This situation can 
be changed by water adsorbates, which may sit either in 
between graphene and the substrate (Fig. [3]c), e) and f) 
or on top of graphene (Fig. [3]b)). In all cases, the elec- 
trostatic dipole moment of the water adsorbates comes 
along with strong local electrostatic fields, which allow 
to shift the impurity bands significantly with respect to 
the graphene bands. As Fig. [3] demonstrates, this shift 
strongly depends on the place of adsorption and the ori- 
entation of the water molecules leading, either to hole 
doping (Fig. [3]c) and f)) or electron doping (Fig. [3]e)). 

This strong effect of water is very general and not lim- 
ited to the examples shown in Fig. [3l Similar effects 
of water can be found for graphene on the (OOOl)-plane 
of Qf-quartz or /3-cristobalite with Q4 defects. Although 
the interaction mechanisms of water, Si02 and graphene 
presented here are not exhaustive, the comparison of wa- 
ter adsorption on perfect free standing graphene to wa- 
ter adsorption on graphene lying on a (defective) Si02 
substrate allows the following conclusions: Perfect free 
standing graphene may have water adsorbates on top but 
its electronic transport properties are insensitive against 
perturbations by the water adsorbates. Single molecules 
will not create any impurity states close to the Dirac 
point. For obtaining doping effects, one needs highly or- 
dered H2O cluster or ice structures. 

The substrate changes the situation completely. The 
dipole moments of H2O adsorbates cause local electro- 
static fields that can shift the substrate's defect states 
with respect to the graphene electrons and cause dop- 
ing. The hybridization of substrate defect states with 
the graphene bands can be reduced by H2O in between 
graphene and the substrate leading to less scattering of 
graphene electrons at defects in the substrate. On the 
other hand, impurity bands can be shifted towards the 



Fermi level by H2O adsorbates, leading to increased elec- 
tron scattering. So, for graphene on a substrate, H2O 
much more likely affects the electronic properties than 
for free standing graphene. The effect of water strongly 
depends on properties of the substrate like the amount 
and type of defects. 

This finding might explain experiments on CNTs 
CNTs on Si02 substrates are much more sensitive to 
water in their environment, than suspended CNTs or 
CNTs separated from the Si02 substrate by poly(methyl 
methacrylate) (PMMA) coating. In experiments simi- 
lar to those in Refs. [1, [lO] these effects can be checked 
for graphene: The creation of electron/hole puddles and 
doping on H2O exposure might be investigated in a study 
with a similar setup as in Ref. [io^. We expect much 
less impact of II2O on the transport properties for free 
hanged graphene. Similarly, a comparison of graphene on 
hydrophobic substrates and hydrophilic substrates like 
PMMA and Si02 in Ref. [3] should yield weak and strong 
response on H2O, respectively. 

These findings can be important for the development 
of graphene based gas sensing devices: Some molecules 
like NO2 lead to acceptor states even without a substrate 
[5], while doping effects due to other molecules like H2O 
strongly depend on the substrate. This can open the 
possibility of selective graphene gas sensors. 
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Figure 3: (Color online) Band structures for graphene on defective Si02 substrates, (a-c) 2x2 and (d-f) 4x4 graphene 
supercells with every second, (a-c), or eighth, (d-f) surface Si atom forming a Q% defect. The corresponding geometries are 
shown in Fig. [T]d-i), respectively. Spin up and down bands are shown at the same time. Contributions at the defect site are 
marked as green fatbands. a) and d) without water adsorbates. b) with water on top of graphene. c), e), and f) with water 
between graphene and the substrate, c), f) H2O dipole moment pointing tilted downwards, e) H2O dipole moment pointing 
upwards. 
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